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Abstract 

The method of squared eigenfunction potentials (SEP) is developed systematically 
to describe and gain new information about Kadomtsev-Petviashvili (KP) hierarchy 
and its reductions. Interrelation to the r-function method is discussed in detail. The 
principal result, which forms the basis of our SEP method, is the proof that any 
eigenfunction of the general KP hierarchy can be represented as a spectral integral 
over the Baker- Akhiezer (BA) wave function with a spectral density expressed in terms 
of SEP. In fact, the spectral representations of the (adjoint) BA functions can, in 
turn, be considered as defining equations for the KP hierarchy. The SEP method is 
subsequently used to show how the reduction of the full KP hierarchy to the constrained 
KP (cKP^m ) hierarchies can be given entirely in terms of linear constraint equations 
on the pertinent r-functions. The concept of SEP turns out to be crucial in providing a 
description of cKP r rn hierarchies in the language of universal Sato Grassmannian and 
finding the non-isospectral Virasoro symmetry generators acting on the underlying 
r-functions. The SEP method is used to write down generalized binary Darboux- 
Backlund transformations for constrained KP hierarchies whose orbits are shown to 
correspond to a new Toda model on a square lattice. As a result, we obtain a series 
of new determinant solutions for the r-functions generalizing the known Wronskian 
(multi-soliton) solutions. Finally, applications to random matrix models in condensed 
matter physics are briefly discussed. 
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1 Introduction 



The primary object of this paper is the Kadomtsev-Petviashvili (KP) integrable hierarchy (for 
comprehensive reviews, see e.g. [|T], [J) and its nontrivial reductions generalizing the familiar 
r-reduction to the SL(r) Korteweg-de Vries (KdV) hierarchy. The KP hierarchy is an infinite- 
dimensional system which admits different alternative formulations and exhibits many types 
of symmetries. Here we are interested in a formulation based on the notion of squared 
eigenfunction potential and the spectral representations of the underlying eigenfunctions it 
gives rise to. Because of its connection to vertex operators, many aspects of this theory 
are algebraic in nature. This allows us to discuss various symmetries of the hierarchy and 
applications to a large class of soliton systems obtained from it via symmetry reduction in a 
systematic manner. 

The KP hierarchy arises as a set of compatibility conditions for the linear spectral problem 
involving the pseudo-differential Lax operator £ and the Baker- Akhiezer (BA) wave function 
ipBAit, A). In recent years, the study of integrable systems of KP type has undergone rapid 
growth due to the applications of the tau-function technique invented by the Kyoto school 
H |], [H . The underlying principle of this method is to represent the relevant soliton potentials 
and Hamiltonian densities in terms of isospectral flows (with evolution parameters (t) = 
(ti = x, £2, • • •)) °f one single function r(i) in such a way that d 2 In r(t) / dt\dt n becomes equal 
to the coefficient in front of D" 1 in the pseudo-differential operator expansion of C n . 

In terms of the r-function, viewed as a function of the infinitely many KP "time" -variables 
(ti = x, t 2 , . . .), the whole KP hierarchy is contained in Hirota's fundamental bilinear identity 
instead of the infinite system of non-linear partial differential equations derived from the 
Sato- Wilson Lax operator approach. The r-function approach bridges the way to several 
physical applications in view of its direct connection to physical objects, such as correlation 
and partition functions. Moreover, it allows a coherent treatment of multi-soliton solutions. 
These solutions of the nonlinear differential equations are generated by the action of the 
Miwa-Jimbo vertex operator X (A,/i) @] (cf. eq. (|3.1| ) below) on the r-function. This vertex 
operator generates an infinitesimal Backlund transformation of the KP hierarchy. The family 
of all vertex operators constitutes a Lie algebra isomorphic to GL (00). The transformation 
r(t) — > (exp(aX (X, n)))r(t) sends a solution of the KP hierarchy into another solution. In 
this way, the action of the infinite-dimensional Lie algebra GL(oo) on the solution space of 
the KP equation is made explicit via the Backlund transformation of "adding one soliton" . 

This powerful formal machinery embeds many other concrete and useful structures rel- 
evant for physical models. Recently a special class of solutions encountered in the matrix 
models of discrete two-dimensional gravity was realized via the imposition of the Virasoro 
type of constraints on the underlying r-function |6|, [7j (see also ||). 

A remarkable feature of the KP hierarchy is the existence of the so called additional 
non-isospectral symmetries which, within the Lax operator formalism, are generated by 
Orlov-Schulman pseudo-differential operators ||. The latter are defined as purely pseudo- 
differential parts of products of powers of the Lax operator C and its "conjugate" A^-operator 
(cf. eqs.( |5.1| ),( |5.2| ) below) and their respective flows form the infinite-dimensional Lie algebra 
W 1+oo 0. In an important recent development Adler, Shiota, and van Moerbeke |13 
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Wi +oc algebra was originally introduced in physics literature |10| as a nontrivial "large N" limit of 
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also [15, |k|) obtained a formula for the KP hierarchy which relates the action of the vertex 
operator X (A, fi) on the r-function to Orlov-Schulman non-isospectral additional symmetry 
flows on the BA wave function. The coefficients in the spectral expansion of X (A, //) act 
as vector fields on the space of r-functions generating Wi +oc algebra as well. Hence, the 
above result relates the Wi +00 algebra acting on r(t) to the centerless Wi +oc algebra of 
non-isospectral symmetry flows acting on the BA function ipBA(t,X). 

There exists an alternative to the r-function method characterization of the KP hierarchy 
evolution equations in terms of (adjoint) eigenfunctions, i.e., functions whose KP multi-time 
flows are governed by an infinite set of purely differential operators {Bk}'^' =l (cf. Def. |2.1| 
below). The latter, by virtue of compatibility of the multi-time flows, have to satisfy the 
so called "zero-curvature" Zakharov-Shabat equations (cf. eq.( |3.3S| ) below). One can then 



show [[H]] that all Bk are obtained as purely differential projections of k-th powers of a single 



pseudo-differential operator C, thus leading to the standard Lax formulation of KP hierarchy. 

Overcoming the formal obstacle of having to define a function via an inverse deriva- 
tive (9" 1 Oevel succeeded in |18| to associate a well-defined (up to a constant) function - 
the squared eigenfunction potential (SEP), to a pair of arbitrary eigenfunction and adjoint 
eigenfunction such that the x-derivative of SEP coincides with the product of the latter 
eigenfunctions. Consequently, a systematic formalism emerged in [I8| for the study of sym- 



metries generated within the KP hierarchy via SEP [[jl|. In a particular example, when 
both eigenfunctions defining the SEP are BA functions, the SEP becomes a generating func- 
tion for the above mentioned additional non-isospectral symmetries of the KP hierarchy 



In the SEP framework, the product of any pair of eigenfunction and adjoint eigenfunction, 
being a x-derivative of SEP, can be viewed as a conserved density within the hierarchy. 
The transition to the important class of constrained KP hierarchies cKP rm [], which are 
Hamiltonian reductions of the general KP hierarchy and whose Lax operators are given in 
eq.( |2.20| ) below, can be effectuated by imposing equality between a linear combination of m 
(m > 1) conserved densities of the above mentioned type and the r-th (r > 1) fundamental 
Hamiltonian density of the KP hierarchy. In such a case, the symmetry generated by SEP 
(called "ghost" flow) is identified with the r-th isospectral flow of the original KP hierarchy. 

The principal merit of our work is to reformulate the eigenfunction formalism of KP 
hierarchy in a new form called squared eigenfunction potential (SEP) method, namely, to 
employ SEP as a basic building block in defining the KP hierarchy. The main ingredient 
of the SEP method is the proof of existence of spectral representation for any eigenfunction 
involving SEP as an integration kernel (spectral density). A link is then provided between the 

the associative, but non-Lie, conformal Wn algebra It turns out to be isomorphic to the (centrally 
extended) algebra of differential operators on the circle jl^l, i.e., the Lie algebra generated by z k (d/dz) n for 
fc£2,«>0. Let us also recall that the "semiclassical" limit (contraction) Wi +oc of Wi +oc is the algebra 
of area-preserving diffeomorphisms on the cylinder [[T(| . 

4 The cKP riTO integrable hierarchies appeared in different disguises from various parallel developments: (a) 
symmetry reductions (2^, [f8|, of the full KP hierarchy; (b) abelianization, i.e., free-field realizations, in 
terms of finite number of fields, of both compatible first and second KP Hamiltonian structures J|4|, |j|; ( c ) a 
method of extracting continuum integrable hierarchies from the generalized Toda-like lattice hierarchies |29] 
underlying (multi-)matrix models; (d) purely algebraic approach in terms of the zero-curvature equations 
for the affine Kac-Moody algebras with non-standard gradations pj . 
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two alternative formulations of the KP hierarchy: one based on the r-function and another 
one based on the SEP method. Furthermore, we apply the SEP method to solve various issues 
in integrable models of KP type and their applications in physics, among them, deriving new 
determinant solutions for the r-function containing the familiar Wronskian (multi-soliton) 
solutions as simple particular cases, as well as identifying them as possible novel types of 
joint distribution functions in random matrix models of condensed matter physics. 

The plan of the paper is as follows. After reviewing the background material in Section 
|2], we first prove in Section ^| that any eigenfunction of the general KP hierarchy can be 
represented as a spectral integral over the BA wave function with a spectral density expressed 
in terms of SEP. When (at least one) of the eigenfunctions is a BA functions, we obtain 
a closed expression for the SEP. When both of the eigenfunctions are BA functions, the 
resulting SEP's are connected straightforwardly to the bilocal vertex operator X (A, \i) acting 
on the r-function. This association leads to a simple alternative proof for the Adler, Shiota, 
and van Moerbeke result |13|, II, |H| [[(J mentioned above. 



A further important observation in Section |3] is that the spectral representation equations 
for the (adjoint) BA functions themselves can be considered as defining equations for the KP 
hierarchy. In other words, our formalism of spectral representations of KP eigenfunctions 
can be viewed as an equivalent alternative characterization of the KP hierarchy parallel to 
Hirota's bilinear identity or Fay's trisecant identity f[3 



Our results in the constrained cKP r m hierarchy follows. In Section f|, using 

the SEP framework we obtain an equivalent description (parallel to the one within the Lax 
pseudo-differential operator approach) of cKP r m hierarchies entirely in terms of r-functions 
only. Namely, we first derive a linear equation for the r-function (eq.( |4.5| ) below), involving 
the bilocal vertex operator X (A, \i) and a set of spectral densities, which uniquely constrains 
the r-function to belong to the cKP rm hierarchy. Furthermore, we provide in Section |] an 
alternative description of cKP r m hierarchies in the language of universal Sato Grassmannian. 

One of the advantages of the SEP approach lies in the fact that it allows for a co- 
herent treatment of the non-isospectral symmetries for KP-type hierarchies. We use this 
feature in Section |5| to demonstrate how the combination of the familiar Orlov-Schulman 
non-isospectral symmetry flows, operating in the full unconstrained KP hierarchy, together 



with certain appropriately chosen additional SEP-generated "ghost" symmetry flows |28], || 
gives rise to the correct non-isospectral Virasoro symmetry generators acting on the space 
of cKP r m r-functionsQ 

The SEP method is applied further in Section [3] to formulate generalized multi-step binary 
Darboux-Backlund (DB) transformation rules of (constrained) KP hierarchies (one-step bi- 
nary DB transformations with SEP have been introduced previously in ref. ^U|). Based on 
these transformation rules and using the fundamental Fay identities, we derive a series of new 
determinant solutions for the r-functions generalizing the known Wronskian (multi-soliton) 
solutions. The binary DB orbits define a discrete symmetry structure for cKP r m hierarchies 
corresponding to a square lattice. We exhibit the equivalence of these binary DB orbits with 
a generalized Toda model on a square lattice. 

Our formalism offers applications to the study of random matrix models in condensed 
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The standard Orlov-Schulman non-isospectral symmetry flows do not preserve the constrained form 



fl2.20|) of cKP r!m hierarchy. 
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matter physics, which we briefly discuss in Section 0, where we also present some discussion 
of the results and an outlook. 



2 Background on General and Constrained KP Hier- 
archies 

The calculus of the pseudo-differential operators (see e.g. H, H) is one of the principal 
approaches to describe the KP hierarchy of integrable nonlinear evolution equations. In 
what follows the operator D is such that [ D , f] = df = df /dx and the generalized Leibniz 
rule holds: 

D n f = jt( n )(d J f)D n -J , ne% (2.1) 

In order to avoid confusion we shall employ the following notations: for any (pseudo-)- 
differential operator A and a function /, the symbol A(f) will indicate application (action) 
of A on /, whereas the symbol Af will denote just operator product of A with the zero-order 
(multiplication) operator /. 

In this approach the main object is the pseudo-differential Lax operator: 

r—2 oo 

L = D r + v J Dj + u i D ^ ( 2 - 2 ) 

j=0 i=l 

The Lax equations of motion: 

^- = [lT,L\ , n = l,2,... (2.3) 

describe isospectral deformations of L. In (|2.3| ) and in what follows the subscripts (±) 
of any pseudo-differential operator A = J2j a jD^ denote its purely differential part (A + = 
J2j>o ajDi) or its purely pseudo-differential part (A_ = J2j>i a -jD~i), respectively. Further, 
(t) = (t\ = x, t%, . . .) collectively denotes the infinite set of evolution parameters (KP "multi- 
time") from (EDI ). 

The Lax operator Q2.2D can be represented equivalently in terms of the so called dressing 
operator W : 

oo 

W = 1 + w n D- n ■ L = WD r W- 1 (2.4) 

n=l 

whereupon the Lax eqs.([2.3|) become equivalent to the so called Wilson-Sato equations: 

^ = L n + /r W - WD n = -L n _ /r W (2.5) 
Further important object is the Baker-Akhiezer (BA) "wave" function defined via: 



^ BA (t,X) = W(e^ x) )=w(t,X)e^ x) ; w(t, A) = 1 + X>^)A~\ (2.6) 
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where 

oo 

£(t,A) = £t„A" ; t x =x (2.7) 

n=l 

Accordingly, there is also an adjoint BA function: 

oo 

r BA (t, A) = W-^e-tM) = w*(t, A)e-^ A ) ; w*(t, A) = 1 + £ (2.8) 



4=1 



The (adjoint) BA functions obey the following linear system: 

L{ip BA (t,X)) = \ r iP B A(t,X) , ^BA(t,\) = L n + /r ^ BA (t,\)) (2.9) 

L*(^ A (t,A)) = \*p BA {t,\) , ^-r BA (t,\) = -(Lr + / r (r BA (t,\)) 

Note that eqs. (|2.3| ) for the KP hierarchy flows can be regarded as compatibility conditions 



for the system (|2.9| ). 

There exists another equivalent and powerful approach to the KP hierarchy based on 
one single function of all evolution parameters - the so called tau-function r(t) || . It is an 
alternative to using the Lax operator and the calculus of the pseudo-differential operators. 
The r-function is related to the BA functions (12. 6f) -(12. 9f) via: 



, W( , A) = «dp>&» = e «« £ E^-Elim A -» (2 , 0) 

^ (t , A) = I(i±ip) e - {< u, = e -«,« g MM A -„ (2 . n) 

T W n=0 T W 



where: 



[A" 1 ] = f A" 1 , ^A- 2 , 1a" 3 , . . .) ; [d] = f A I * 1 « . . ( 2 .12) 
L J V 2 ' 3 ' y 1 J {dtx' 2 dt 2 3 dt 3 : J K 1 

and the Schur polynomials p„(t) are defined through: 

oo 

e z]i>i A '*i = ^ X n p n (ti,t 2 , . . .) (2.13) 

n=0 

Taking into account ( 2.1 0| ) and ( |2.6| ), the expansion for the dressing operator (|2.4j) becomes: 



w = j^ Pn( [d])r(t) D _ n ^ t . e . Wl ( t ) = ResW = -d x \nr(t) (2.14) 

The (adjoint) BA functions enter the fundamental Hirota bilinear identity: 

d\ij BA (t,X)r BA (t',X) = (2.15) 
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which generates the entire KP hierarchy. Here and in what follows integrals over spectral 



A=0- 



parameters are understood as: / dX = § = Res 

Let us also recall that the KP hierarchy possesses an infinite set of commuting integrals 
of motion w.r.t. the compatible first and second fundamental Poisson-bracket structures 
whose densities hi-\ = jResL'/ r are expressed in terms of the r-function ( |2.10| ) as: 



d x ^-]nr(t) = ResL l/r (2.16) 
oti 

Below we shall be particularly interested in reductions of the full (unconstrained) KP 
hierarchy (|2.2p . In this respect, it turns out that a crucial role is played by the notions of 
(adjoint) eigenfunctions of KP hierarchy. 

Definition 2.1 The function $ (ty) is called ( adjoint) eiqenfunction of the Lax operator L 
satisfying Sato's flow equation t \2.3j ) if its flows are given by the expressions: 



£ = Lf($) ; | = -(Lf r W (2-17) 



for the infinitely many times t^. 



Of course, according to fl2.9|) the (adjoint) BA functions are particular examples of (adjoint) 
eigenfunctions which, however, satisfy in addition also the corresponding spectral equations. 

In what follows a very important role will be played by the notion of the so called 
squared eigenfunction potential (SEP). As shown by Oevel fIBfl , for an arbitrary pair of 
(adjoint) eigenfunctions $(t), *ff(t) there exists the function S ($(t), ^(t)), called SEP, which 
possesses the following characteristics: 

d S ($(t), #(*)) = Res (D-^iL^+QD- 1 ) (2.18) 



dt,_ 

The argument in |Jl8[] , proving the existence of S* ($(£), \I/(t)), was built on compatibility 



between isospectral flows as defined in eq. ( |2.18| ) and (|2.17|) 



In particular, for n — 1 eq. (|2.18|) implies that the space derivative (recall d x = d/dti) of 
S (&(t), ^{t}) is equal to the product of the underlying eigenfunctions: 

d x S = (2.19) 

Remark. Eq.( |2.18D determines S ($(t), ^(t)) up to a shift by a trivial constant. 

From eqs.( |2.18] )- (|2.19|) one sees that $(t) ^(t) is a conserved density of the KP hierarchy. 
This fact has a special significance for the reduction of the general KP hierarchy to the 
constrained cKP rm models (see below). 

Definition 2.2 The constrained KP hierarchy (denoted as cKP rimy ) is given by a Lax operator 
of the following special form: 

r—2 m 

L = L r , m = D r + Y, u t D l + £ ^ZT 1 * a (2.20) 

1=0 o=l 
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One can easily check that the functions $ a , \l/ a , entering the purely pseudo-differential part 
of L r>m ( |2.20|) , are (adjoint) eigenfunctions of L r m according to Def.|2.1|. 

The purely pseudo-differential part of arbitrary power of the cKP r r7l Lax operator ( |2.20| ) 
has the following explicit form 



m k—1 



L l = E E L^-^D- 1 (L*Y (tf a ) (2.21) 
0=1 3=0 

Formula ( p.21|) can easily be derived from the simple technical identity involving a product 
of two pseudo-differential operators of the form fiD~ x gi , i = 1, 2 : 

f 1 D- 1 g 1 f 2 D- 1 g 2 = frS^, gi)D~ x g 2 - hD^Sfo, 9l )g 2 (2.22) 

where /j, g^ are pairs of (adjoint) eigenfunctions of some KP Lax operator, with S(-, •) being 
the corresponding SEP. 

Note, that in agreement with eq.( 2.22|) the expression L(<$> a ) in ( 2.21|) with L as in ( |2.20|) 
explicitly reads: L($ a ) = L + ($ a ) + J2l=i &bS (<& > ^b)- This definition extends naturally to 
higher powers of L acting on $ a as well as higher powers of L* acting on \l/ a . Moreover, one 
can easily show that L l ($ a ) and L* k (^? a ) are (adjoint) eigenfunctions of L fl2.20|) as well. 

For three pseudo-differential operators X t = fiD~ x gi ,i = 1,2,3 the associativity law 
(XiX 2 ) X 3 = X\ (X 2 X 3 ) implies the following technical Lemma: 

Lemma 2.1 The squared eigenf unction potential S(-, ■) satisfies: 

S(f, g)S(h, k) = S (h, kS(f, g)) + S (fS(h, k),g) (2.23) 
for (adjoint) eigenfunctions f, g, h, k. 



3 Spectral Representation of KP Eigenfunctions 

Consider the bilocal vertex operator |4j : 

X (A, n) = - : e*» : : e~ S ^ := L^-M^-eM) K^"^ 
A A 

where: 
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^.rt-flt-rL^L, 1^ )^i + S(\,fi) (3.1) 



W^-E^ + Ey^V (3-2) 
1=1 1=1 1 at l 



£ (t, A) is as in ( |2.7|) , the columns : ... : indicate Wick normal ordering w.r.t. the cre- 



ation/annihilation "modes" t\ and -M-, respectively, and the delta-function is defined as: 



, . 1 1 1 1 . . 

SM = (3.3) 
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An equivalent representation for X (A,/i), using Wick theorem, reads: 



X(\, fi) = — — : e § W~ § M ;= _J_ e «*,M)-€(*,A) J2T H*-'-»- l )£i for ^1 < | A | (3.4) 

A — LI A — yU 

The vertex operator X (A, //) can be expanded as follows: 

* (A, /i) = ^- £ f A— ^j^rW^ (3.5) 

where the operators span Wi +00 algebra. 

From the standard representation for the (adjoint) Baker- Akhiezer wave function (|2.10|) , (|2.11|) 
in terms of the r-function we deduce the identity: 

= —if> BA (t, t j l )if>* BA (t-\ii- 1 ],\)+6{\,iJ l ) (3.7) 

Now, recall the Fay identity |13| for r-functions: 

(s - Sl )(s 2 - s 3 )r{t + [s ] + [ Sl ])r(t + [s 2 ] + [s 3 ]) + cyclic(l, 2, 3) = (3.8) 

which, in fact, is equivalent to Hirota bilinear identity (|2.15 ). In what follows, we shall often 
make use of a particular form of (|3.8|) upon setting Sq = 0, dividing by SiS 2 s 3 and shifting 
the KP times (t) -> (t - [s 2 ] - [s 3 ]) : 

(sj 1 - s, 1 ) r(t+ [ Sl ] - [s 2 ] - [s 3 ]) r (t) + (sf 1 - s 2 l ) r{t- [s 2 ]) r(t+ [s,] - [s 3 ]) 

+ («3 1 - ^r 1 ) r(t- [s 3 ]) r(t+ [ 8l ] - [s 2 ]) = (3.9) 

Especially, making identification si = s 2 = z~ x and S3 = A" 1 in (|3.9| ) and using ( |2.10| )- 
( P-llj ), we arrive at the following useful Lemma: 



Lemma 3.1 The truncated Fay identity is equivalent to the following bilinear identity 
for (adjoint) BA functions: 

^£,(V^,A)^*- [A" 1 ],//)) = -- z i)BA{t,XW BA (t-[z-\v) (3-10) 

where A z is the shift- difference operator acting on functions depending on the variables t = 
(ti,t 2 , ...) as follows: 

= e£r - l , A z f(t) = f (t - [z- 1 ]) - fit) (3.11) 
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The Fay identity ( j3.8|) is also known in its differential version: 



d ( r (t + [A' 1 ] - [^}) \ =(x _ u) ( r (t + [A' 1 ] - [^]) _ r (t + [A' 1 ]) r (t - ^]) \ 
X \ r(t) J { N [ r(t) r(t) r(t) ) 

(3.12) 

Using ( |3.4|) and multiplying both sides of (|3.12|) by exp { — £ (i, A) + £ (t, //)} we can rewrite 
it as: 

AS. ( Ar(A ;^ T(t) ) = ~1>BA(t, Wba (t, A*) (3-13) 
or, equivalently, using (|3.6|) and 



3s f-^Bxft A) VfiA (t + [X Xfifj = ^Bx(*> A)^BA (*, /*) 
^f^BA^//)^^-^" 1 ]^)^ = r B A(t,WBA(t,v) (3.14) 

Let \I/ be a pair of an eigenfunction and an adjoint eigenfunction of the general KP 
hierarchy. Our main statement in this section is: 

Proposition 3.1 Any (adjoint) eigenfunction of the general KP hierarchy possesses a spec- 
tral representation of the form: 



$(t) = J dX ip(X) ip BA (t, A) ; (f) = J dX <p*(X) r BA (t, A) (3.15) 
with spectral densities given by: 

¥>(A) = \Pba A) $ (f + [A" 1 ]) ; ^(A) = ^A(t',A)vI/(t'-[A- 1 ]) (3.16) 
where the multi-time t' = (t^t^, • • •) taken at some arbitrary fixed value. In other words: 

$(t) = Jd\iJ> BA (t,\)^ BA (lf,\)$(lf+[\- 1 ]) (3.17) 
9(t) = J dXr B A(t,X)\^BA(t',X)^ (f -[X- 1 ]) (3.18) 

are valid for arbitrary KP (adjoint) eigenfunctions \l/ and for an arbitrary fixed multi-time 
t' . Furthermore, the r.h.s. of ( \3.11\ ) and (\3.18Q do not depend on t' . 

We will proceed proving the above proposition in two steps. First, let us assume that the 
(adjoint-)eigenfunctions indeed possess a spectral representation of the form fl3.15p with some 
spectral densities ip(*'(X) . In such case the statement of the proposition is contained in a 
much simpler Lemma: 



Lemma 3.2 For (adjoint) eigenfunctions possessing the spectral representation ( S.lfy) their 
respective spectral densities are given by ( \3.16j ). Consequently, in this case the equations 
( $.1T[ ) and ( $.1$ ) are valid too. 
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Proof. Using the spectral representation ( j3.15| ) for $ (t f + [A x ]) and substituting it into 
the right hand side of (|3.17|) , we get: 



dX J dyL <p(fjl) if> BA (t, \)^P*ba (t', A) IpBA (t + [A" 1 ], n) (3.19) 
Recalling (|3.7p we can rewrite ( |3.19|) as: 

dX J dfnp(fi) ip B A(t,\) 1—ipBA {t',fj)^ B A (f - \m~ X ]^) +S(X,fl)j 



dX<p(X) ^ BA (t,X) = 



(3.20) 



where use was made of the fundamental Hirota bilinear identity (|2.15|) . The t'-independence 
of the r.h.s. of ( |3.17| ) and ( |3.18| ) will be demonstrated explicitly in the course of proof of 
Prop.|3j] given below. □ 

We are now ready to take a final step of the proof of Prop. |3.l| and extend the result of 
Lemma (|3.2|) to arbitrary KP (adjoint-)eigenfunctions without assuming existence of a spec- 
tral representation (|3.15|) . To this end we need to recall the notion of SEP (|2.18|) -( f2.19|) . 

Let S ($(t), ^ba^i A)) be the SEP associated with a pair of eigenfunctions <£>(£) and 
*p* BA {t,X), i.e. d x S($(t),il>* BA (t,\)) = ${t)ij)* BA (t,\). Define now: 



dXip BA (t, A) S ($(t'),ip BA (t', A)) 



(3.21) 



We first observe that .9$ (t, t') /dt' n = due to eqs. (glSD and fgTT5|) . Hence $ (t, f) = $ (t) 
does not dependent on the multi-time t! . Moreover, it is obvious from the definition ( |3.21| ) 
that <3> (t) is an eigenfunction possessing spectral representation of the form ( |3.15|) and, 
therefore, satisfying the conditions of Lemma (|3.2|). Consequently, according to ( |3.17| ), we 
have: 

S(f) = / dX^ BA (t,X)jr B A(t',X)^ (f + [A- 1 ]) (3.22) 



Agreement between ( |3.21| ) and ( |3.22| ) requires that their respective A-integrands may differ 
by at most a term proportional to the Hirota A-integrand in eq. ([2.15| ). The latter implies 
the fulfillment of the following identity: 



jf BA (t, A)$ (t + [A- 1 ]) = -S ($(t), r BA (t, A)) + Ar BA (t, A) 



(3.23) 



where A is some differential operator w.r.t. t = (ti,t 2 , . . .). Acting with d x on both sides of 
( |3. 23| ) and using identity (p . 1 4|) together with ( |3.22|) , we conclude that A has to satisfy 



$ (t) - $ (t) 



d x Ar BA (t, A) 

r BA (t, a) 



(3.24) 



The left hand side of ( |3.24|) is a KP eigenfunction, independent of the spectral parameter A. 
This forces A = and consequently $(£)=<& (t). This concludes the proof of Prop.|3.1|. □ 
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Corollary 3.1 Taking into account Prop. 3.1, eqs.( 3.13] )-( 3.14 ) imply the following rela- 
tions: 



S (V>Ba(*, A), V>BA 

S($(t),^ A (t,A)) 
S{i/> BA (t,\),*{t)) 



-\*ba (t + [A" 1 ], n) r BA (t, A) = - X{X £l T{t) (3-25) 



-l/>BA(t,\)ll>*BA f*-^- 1 ],^' 



r(t) 
*(//,A)r(t) 
r(t) 



--TP* BA (t,\)<f>(t+[\- r 



8(/i, A) 
(3.26) 
(3.27) 

(3.28) 



where $, \& are arbitrary (adjoint- jeigenj 'unctions and £(•,•) ore i/ie corresponding squared 
eigenfunction potentials. Moreover, we also have the following double spectral density repre- 
sentation for the SEP S ($(*), #(£)); 

S ($(*),¥(*)) = -// dAd//^(A)^(/x)^ A (t, A)^(t + [A" 1 ],//) (3.29) 



Taking into account ( |3.27| )-( p.28| ), the spectral representations ( |3.17| )-( |3.18| ) become: 

m = -Jd\iJ BA (t,\)sm'),r BA (t'^)) (3.30) 

*(*) = fd\p BA {t,\)S(il> BA (lf,\)Mt)) (3-31) 
Remark. Note that the expressions ( 3.30 )-( 3.31|) applied for (adjoint) BA functions yield: 

1pBA(t,X) = ~ J dflll>BA(t,ll)S (V , flA(*',A),V'flA(* , >^)) 

^(*>A) = fdfAil>% A (t,riStyBAW,tt),i>%AV,X)) (3-32) 



which shows that the SEP 5 (ipBA(t', A), i[> BA (t', /i)) can be identified with the Cauchy kernel 
for each fixed KP multi-time t' (cf. also ]2U| and references therein, where the above SEP 



was previously introduced in the context of Riemann factorization problem, as well as |21 
for related discussion within the dispersionless KP hierarchy). 



Remark. Going back to the spectral representation eqs. (ft.l7|) - (|3.18 ), valid for any eigen- 
function of the general KP hierarchy, we observe that they can be rewritten as evolution 
equations w.r.t. the KP multi-time of the following form: 



, U(t,t') = |rfA^(t,A)^ A (t',A)eSr^- i ^; (3.33) 
U*(t,t') = JdXr B A(t,^BA(t',X)e-^T^-'9/dt' l ( 3.34) 
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One can readily verify that: 



U(t, t) = l , ET 1 (t, t') = U(t', t) , U(t, t') = U(t, t") U{t", t') (3.35) 

±U(t,t') = LfU(t,t') , ^U(t,t') = -Lfu(t,t') (3.36) 
oti ot\ 



From ( p.35| )- p.36| ) we deduce that the evolution operator U (t, t') (|3.33|) can be formally 
written ClS db path-ordered exponential: 

{OO „\ Jj. 
J2J ds-^L l i r (t(s))j ; t k (0)=t' k , t k (l)=t k , k = l,2,... 

(3.37) 

which precisely agrees with the formal solution of the differential evolution eqs.( |2.17| ) for 
the KP eigenfunctions. The r.h.s. of ( [3.37D is independent of the particular path {t k (s)} 



connecting the points t! and t in the space of KP multi-times due to the "zero-curvature' : 
Zakharov-Shabat equations: 



d l/ r d u/r 



jk/r jljr 



(3.38) 



Thus, our SEP method allowed us to find the explicit expression (r.h.s. of the second 
eq. (|3.33|) ) for the formal path-ordered exponential ( |3.37|) . 

Now, it is worthwhile to observe that we can revert the logic of our procedure above, i.e., 
instead of starting with Hirota bilinear identity ( [2.15| ) (or, equivalently, with Fay identity 
( |3.8D ) as defining the KP hierarchy and deriving from them the spectral representation 
formalism ( |3.17| )-( |3.18| ) (or ( |3.30| )- (|3.31| )) for KP eigenfunctions, we can take the spectral 



representation eqs. ( p.30| ) (|3.31| ) as the basic equations defining the KP hierarchy. Namely, 
we have the following simple: 

Proposition 3.2 Consider a pair of functions tp(t, \),ip*(t, A) of the multi-time (ti,t2, ■ ■ 
and the spectral parameter A of the form ip^(t,X) = e^*'^ J2JLo Wj (t)X~^ with Wq*^ = 
and ^(t, A) as in ( \2. Let us assume that ip^*\t, A) obey the following spectral identities: 



4j(t,X) = - dfi^{t,n)S{t';X,fi) , ip*(t,X) = dni>*(t,n)S(1?',ti,\) (3.39) 



for two arbitrary multi-times t and t' , where by definition the function S(t;X,fi) is such 
that ^-S(t; A, /i) = ip(t, A) ip*(t, /i). Then, eqs. ( \3. 39\ ) are equivalent to Hirota bilinear iden- 
tity ( $.13i) and, accordingly, ip^*\t,X) become (adjoint) BA functions of the associated KP 
hierarchy. 



To see that eqs. (|3.39|) imply Hirota identity ( |2.15|) , it is enough to differentiate both sides of 



( fOg ) w.r.t. t[: = dip(t, X)/dt[ = —ip(t', A) / dpLip(t, fx) ip*(t', fx). The proof of the inverse 
statement of the equivalence, namely, that Hirota bilinear identity (|2.15| ) imply the spectral 
representation eqs. (|3.39|) , is contained in the proof of Prop.|3~T| above. 
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Using (EnqHEmP , eqs.(gg7|)-( ggg) can be rewritten as: 

Ar(t) 

rCt-^D^ft-tA- 1 ]) e(t>A) 
Ar(t) 



-5($(t),^(t,A)) 
S'(V'fl J 4(t J A),*(t)) 



(3.40) 
(3.41) 



Remark. Spectral representations for eigenfunctions (|3.30|) -( |3.31| ) as well as identities 
( p.40| )- (|3.41 ) were obtained in a similar form in [31] for the particular case of constrained 
cKP r m hierarchies. Let us specifically emphasize, that all main equations of the present SEP 
method (l3l5D - (gl8D , (T3^ - flOlD and (ggjjH ggg), derived above, are valid within the 
general unconstrained KP hierarchy. 

Acting with space derivative d x on both sides of (|3.40 )-( 3.41 ) and shifting the KP time 
arguments, we get: 

^~ [A " ]) - 1 + A-*ln*(t) = A-^ln r ^/f^ (3.42) 



«&(*) 



T(t) 



*(*) r(f) 



(3.43) 



which were obtained in || by studying the way the r-function transforms under Darboux- 
Backlund transformations. Taking into consideration that: 



$(t) A«"!$(t) 



$(t) 



(3.44) 



with Pn(-) being the Schur polynomials (|2.13| ), we find that eq.( |3.42| ) is a generating equation 
for the following set of equations upon expanding in powers of A" 1 : 



Pn[ 



$(*) = v n {t)$(t) ; n>2 
v n (t) = p n ^(-[d})d\nr(t) 



(3.45) 



Note that v n (t) are coefficients in the A-expansion of the generating function v(t, A) |32| : 

oo 

v(t, A) = Y J v n+l \ n = d x In ipBA{t, A) - A = A x d x \nr(t) (3.46) 



71=1 



where in obtaining the last equality we again used eqs.( 2.1d| )— ( 2.1 1| ) and notation ( |3.11| ). We 
will later need a slight generalization of ( |3.46|) : 



,(') 











(t,X) = J2^(t)X n ^—\mlj BA (t,X)-X l = A x —\nr(t) ; 1>1 



n=l 



dti 



(3.47) 



Clearly cr$(t) = p n (— [d]) d/dt t lnr and v n (t) = cr^l^t) , n > 2. The coefficients a$ enter 
the basic identity for the KP Lax operator fl2.2|) : 



(3 At 



CXJ 



n=l 
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Remark. Eqs. (|3.45| ) are, clearly, valid for an arbitrary eigenfunction $ of the full KP 
hierarchy On the other hand, in ref. (see also [p3fl ) eqs. (|3.45|) were presented for the 
special case of $ = ipBA{t, A) as relations equivalent to the standard KP evolution equations 
d^BA(t, A)/<9t n = {L n / r ) + il)BA{t, A). In fact, as shown in [fffil , plugging the BA wave function 
= ipBA{t,fi) into eq.( |3.42| ) one easily recovers the differential Fay identity ( |3.12j ). 

E23l M . Let d n be a 



We now define the "ghost" symmetry flows generated by the SEP [jn|, 
vector field, whose action on the Lax operator L and, accordingly, on the dressing operator 
W, is induced by a set of (adjoint) eigenfunctions <3> a , \l/ a , a G {a} through: 



d a L= [ ]T ^LT 1 ^,! 

aS{a} 



J2 ®aD-^ a )W 

ae{o} 



(3.49) 



As shown in |18| , the corresponding action of the above "ghost" flows on the (adjoint) 
eigenfunctions 



a£{a} 



d a v= J2 s($ a ,y)y a 

a&{a} 



(3.50) 



is compatible with the isospectral evolutions of $, Furthermore, it is easy to see that 

d a S = J2 S($,^ a )S($ a ,^) (3.51) 

ae{a} 

is compatible with eq.( |3.50| ) 

If d/3 W = (j2be{p} ^bD^^b) W defines some other "ghost" flow and both flows d a and 
dp satisfy ( |3.50| ), then: 

[d a ,d p ]W = (3.52) 
as follows from the technical identity (|2.22|) . Equations ( |3.50|) and ( |3.52 ) can be compactly 



raa 



expressed by an identity [d a - EaeM ®aD l W a , dp - J2be{/3} ®bD ^ 

Define now Y(X, /i) = ipBA (t, /i) D~ 1; ip BA (t, A) (cf. ref.[[L5|]) to be pseudo-differential oper 
ator inducing a ghost-flow d(\ tfl )W = Y(\,fi)W according to (|3.49|) . In this case the "SEP' 
symmetry flow is generated by an infinite combination of W 1+oc algebra generators 
Then, according to eq.( |3.50|) the action of this flow on the BA wave function is given by: 



Y(X, ll)(lp BA (t, z)) = d( X)ll )(lf)BA(t, «)) = ^BA (t, /J.) S (lp B A(t, z), 1p* BA (t, A)) 



(3.53) 



Further, let us also define the action of the vertex operator X (A, fi) on the BA function 
ipBA(t,z) as generated by its action (as a vector field) on the ratio of r-functions entering 

* (A, p) *Pba (t, z) = e«'." T(f) ^ (A ' "> T(t ~ '-•"'»- f - |2 " I)?(A - "> T(<) (3.54) 

r (t) 

The latter, upon using the shift-difference operator ( 3.11|) , can be written as: 



X (A, (j) 1pBA (t, z) = IpBA (t, z) A z 



T{t) 



(3.55) 
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Let us stress that, according to ( |3.53| )- (|3.55| ), Y(X, /i) acts on the BA function as a standard 
pseudo-differential operator, whereas X (A, \i) acts on it as a shift- difference operator. 

Now, the above results allow us to give a simple straightforward proof of the following 



version of the Adler-Shiota-van-Moerbeke proposition |TJ, [L5| . It provides the connection 
between the form of the non-iso 'spectral ("additional" ) symmetries of KP hierarchies acting 
on the Lax operators and BA functions 0, on one hand, and their respective form when 
acting on KP r-functions, on the other hand. 

Corollary 3.2 With definitions (\3. 53\) and ( \3. 54\ ) it holds: 



X (A, n) ip B A(t, z) = Y(X, n){ipBA{t, z)) 



(3.56) 



Proof. Indeed, applying ( ^.7| ) and Lemma |0| to the r.h.s. of Q3.55| ), the latter equation 
can be rewritten as: 

X (A, fi) IpBA (t, Z) = -IpBA (t, Z) IpBA (t, (J.) ^* BA {t ~ [z~\ A) = Y(X, fi) (ip B A (t, z)) (3.57) 

z 

where in order to arrive at the last equality use was made of ( |3.26| ). □ 

In the literature one often comes across the vertex operator defined as X(X, fi) 
: exp (§(X) — 9(fi)J : = (A — fi)X (A,//). In such a notation the expression ( |3.56j ) becomes 

X(X, fi) = (A - fi)Y{X, n) as in 0, |S|. 

We conclude this section by proving the following important property of SEP: 



Lemma 3.3 Under shift of the KP times, the squared eigenfunction potential obeys: 

1 



s ($(t - [a- 1 ]), *{t- [a- 1 ])) - s mi *(«)) = -j*wa - [a- 1 ]) 

S ($(t + [A- 1 ]), *(t + [A- 1 ])) - S ($(t), = -U(t + [A" 1 ])^) 



(3.58) 
(3.59) 



Proof According to ( ggg ) and ( ggg ) : 

A 2 S($(t),^(t)) = / / dA d// 0*(A)0(/i) (Vba(*,A),Vba (*^)) 



while from eq. (|3.10| ) we find that: 



A 2 S (^ba(*, A),^ A (t,//)) = — ^ba(*, A)^ A (t - [z \,fi) 



Inserting the above identity back in (|3.60| ) gives ([3.58 ) . □ 

After expanding identities (|3.58|) and (|3.59|) in power series w.r.t. A we obtain: 



Ps (-[d})s($(t)Mt)) = -mPs-i(-[9]Mt) 

Ps ([d])s^(t)^(t)) = mps-immt) , 

where p s (-) are the standard Schur polynomials (|2.13|) . 



8 = 1,2, 



(3.60) 



(3.61) 



(3.62) 
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4 Constraints on cKP rm Tau- Functions. Grassmannian 



Interpretation 



From now on we concentrate on studying the class of constrained cKP r;m hierarchies for 
which we have: 



(4.1) 



a=l 



according to eq. (|2.20| ). We first note that the cKP rjm BA function satisfies, according to 
Ql.lD , the following spectral equation: 



L r , m i>BA(t, A) = X r ipBA (t, A) = (L r , m ) + ^ BA (t, A) + £ $ a (t) S (^ A (t, A), tf a (t)) (4.2) 



a=l 



Due to eq. (|3.28 ), the latter can be cast in the following form: 



X r ^BA(t, A) = (L r , m )+^ BA (t, A) + £ y*a(t)^a (t - [A" 1 ]) fo A (t, A) 

a=l ^ 

m 

= (L r , m )+^(*, A) - [S ($«(* - [A" 1 ]), (i - [A- 1 ])) - 5 ($ a (t), a (*))] ^(*, A) (4.3) 

a=l 

where the second equality in (fOf) follows from (|3.58|) . Recalling relation ( p. 47] ) we find that 
dr(t)/dt r = Y%=x S ($„(*), V a (t)) r{t). Similarly, using the spectral identity L n rr J) BA {t, A) = 
X rn ip BA (t, A) and taking into account relation (|2.21|) we obtain the following set of differential 
equations for the cKP r rn r-function: 



d 



8t r 



■r(t) 



m n—1 



EE^( i,wl i $ «). £,i (*.: 



.a=l i=0 



it) 



(4.4) 



Using the differential Fay identity ( |3.12 ), eqs.( f4.4|) can be equivalently written in the form: 







Bt 



_d_ 



„ „ m 

dXd^ (X r - //y 1 £ <Pl(X)<Pa(fi)X (A, /i) 

0=1 



r(t) = 



(4.5) 



where y a *^(A) are the "spectral densities" of the (adjoint) eigenfunctions & a (t), ^ a (t) entering 
the pseudo-differential part of the cKP rm Lax operator fl2.20|) , and also we have used the 
identity: 

X(\,n) =(fi l -X l )x(X,fi) (4.6) 

Thus we arrive at the following statement providing an alternative definition of cKP rm hi- 
erarchies intrinsically in terms of r-functions: 

Proposition 4.1 Reduction of the full KP hierarchy (j2.2\) to the cKP r m hierarchy in terms 
of Lax operators ( 2.2Q) is equivalent to imposing eqs.( \4-fy) as constraints on the pertinent 
t -functions, where (p^*\X) are "spectral densities" of KP (adjoint) eigenfunctions given as 
in eqs.( \3.16j ). 
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Let us now translate eg. ( f4.5|) into the language of universal Sato Grassmannian Qr || j33| . 
Consider the hyperplane W G Qr defined through a linear basis of Laurent series {/^(A)} in 
A in terms of the BA function as generating function F(t, A) : 



W = span</i(A) ,/ 2 (A),...) 



A(A) 



dx 



x=t 2 =t 3 = 



F(t, A) = ipBA(t, A) 



(4.7) 



In case of the standard r-th KdV reduction, where the corresponding Lax operator C = 
D + YLT UiD~ l satisfies C = C r + , the latter constraint translates to the Grassmannian 
language as A r W C W §. 

Our aim now is to express the cKP rm constraint (|4.1|) (cf. fl2.20|) ) in the Grassmannian 
setting. We find from ( |4.3|) that the generating function F'(t, A) : 



Pn {-[d])s^ a {t)^ a {t)) 



a=l n=l 



A" 



i[> BA (t, A) = {L r , m ) + ^ BA {t, A) (4.8) 



defines via ( |4.7|) a point W of Sato Grassmannian Qr : 

W = span(F'(0, A) , d x F'(0, A), d 2 x F'(0, A), . . .) (4.9) 

which coincides, because of the second equality in (|4.8j), with the original point W defined 
through F(t, A) = ip B A{t, A) (p7|) . Thus, we haveQ : 

Proposition 4.2 Let S (<3> a (t), ^ a {t)) , a = 1, . . . ,m , be m squared eigenfunction potentials 
( \2.1$ ) where $ a , \I/ a are (adjoint-)eig en functions of the general KP hierarchy ( \2.Q ). Then, the 
reduction of (2A) to the cKP r m hierarchy l \2.2(] ) can be equivalently expressed as a restriction 
of Qr to a subset whose points (hyperplanes) W &4-V are subject to the following constraint: 



X r + J2 A A S($ a (t),tt a (t)) 



a=l 



W c W 



(4.10) 



with A\ as in ^3.11 ) and S ($ a (i), \l/ a (t)) being given by ( 3.2L ) in terms of the generating 
function fiflj ) ofW. 



5 Non-Isospectral Virasoro Symmetry for cKP r i r-func- 
tions 

The conventional formulation of additional non-isospectral symmetries for the full KP inte- 
grable hierarchy || [Ll| is not compatible with the reduction of the latter to the important 



class of constrained cKP rim integrable models. In refs. p8|, || we solved explicitly the problem 
of compatibility of the Virasoro part of non-isospectral symmetries with the underlying con- 
straints of cKP rjm hierarchies within the pseudo-differential Lax operator framework. Our 

6 For a different criteria characterizing cKP rjm hierarchies within the Sato Grassmannian framework, see 
refs.il M. 
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construction in [28, g] involves an appropriate modification of the standard non-isospectral 
symmetry flows, acting on the space of cKP rjm Lax operators, by adding a set of additional 
"ghost symmetry" flows (of the type appearing in eq.( |3.49| )). In this section, we derive the 
explicit form of the action of the correct modified Virasoro non-isospectral symmetries as 
flows on the space of cKP rm r-functions. Note that the corresponding result for the full 



unconstrained KP hierarchy has been previously obtained in [14, 15, 16 



To this end, let us first recall that the standard additional (non-isospectral) symmetries [Q, 
16 1 are defined as vector fields on the space of general KP Lax operators (|2.2|) or, alternatively, 
on the dressing operators (|2.4D, through their flows as follows: 



dk, n L 



(M n L l 



[M n L h 



+ nM n -"L 



n—l t k 



d k , n W 



(M n L l 



Here M is a pseudo-differential operator "canonically conjugated" to L such that: 



L , M 



\ 



d_ 

dti 



M 



-l/r 



M 



W 
(5.1) 

(5.2) 



Within the Sato- Wilson dressing operator formalism, the M-operator can be expressed in 
terms of dressing of the "bare" operator: 



M (0) = £ ^ fiD l-r = X{ 



1>1 



„•) • Y.——'' i ,)l 
i>i 1 



x, 



i 



(r) 



(5.3) 



i=i 



conjugated to the "bare" Lax operator = D' 



The additional symmetry flows ( |5 . 1| ) commute with the usual KP hierarchy isospectral 
flows given in (|2.3| ). However, they do not commute among themselves, instead they form a 
centerless Wi +oc algebra (see e.g. |16|). One finds that the Lie algebra of operators Bk, n is 
isomorphic to the Lie algebra generated by —z k (d/dz) n . Especially for n = 1 this becomes 
an isomorphism to the centerless Virasoro algebra dk,\ ~ —Ck-i, with [Ci, Ck] = (l — k)Ci+k- 

As demonstrated in |28|, || , the conventional non-isospectral flows (|5.1| ) do not preserve 
the space of cKP rm Lax operators given by fl2.20| ). In particular, for the Virasoro non- 
isospectral symmetry algebra the transformed Lax operator Bk,iL belongs to a different class 



of constrained KP hierarchies cKP 



r,m(k— 1) 



(when k > 3). The solution to this problem is 



provided by the following [28 



Proposition 5.1 The correct non-isospectral symmetry flows for the cKP r 
( 2.2Q) , spanning the Virasoro algebra, are given by: 



(ML k ) 



_(_ vi 1 ) 



hierarchies 



(5.4) 



i.e., with the following isomorphism £k-i ~ — [ML } 
symmetry generating operators (cf. (\3.4Q) ) defined as: 



+ X« 



/.-i. where x l-i 



(1) 



m k— 1 

EE i / 

t=l j=0 



(A; - 1) ) L h ~ l ~ j (®i)D~ x (L 



k > 1 



are ghost- 



(5.5) 
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Since (auto-)Darboux-Backlund transformations of cKP r m hierarchies (see next section) 
play a fundamental role for finding exact solutions, as well as in establishing the link be- 
tween cKP r m integrable models and (multi-)matrix models, it is natural to impose the ad- 
ditional condition of commutativity of the non-isospectral symmetries with the Darboux- 
Backlund transformations. The latter condition was shown in refs.p8|, |J to be satisfied only 
by the subclass cKP r l of constrained KP hierarchies (it is precisely cKP r l hierarchies which 
provide the integrability structure of discrete multi-matrix models M). Therefore, in the 
rest of this section we restrict our attention to cKP r l models. 

Consider the modified non-isospectral Virasoro symmetry flows 
ing operator of cKP r i hierarchy: 



d* k W = — [ML k 



w + x { k \w 



acting on the dress- 
(5.6) 



Taking the operator residuum on both parts of (|5.6|) we obtain: 



= ^KLifV + 



2r 



■fc-2 

E 

■3=0 



(fc-2) -j ) S(L k - 2 ^^),L* j ^) 



r(t) 



(5.7) 



In deriving eq.(|5.7|) we used the expression for X^_ x 
identity (|3.12|) as well as: 



together with the differential Fay 



Res [M l L k 



Res) 



(M l L k ) (ij BA (t,X))r B A(t,X) 



(5i 



-cL 



T{t) 



Res) 



r(t) 



d x 



T 



In the chain of the identities in (|5.8|) we took into account Dickey's formula for [M l L k 

fT5| (first equality in (|5~8|)), eq.( |3.13| ) (third equality in (|5.8|)), and formula ( |3.5| ) for X (A,/i) 
to arrive at the last equality above. The Virasoro operator in the first term on the r.h.s. 
of ( |5.7| ) comes from the standard Orlov-Schulman non-isospectral symmetry flow and reads 
explicitly (for k > —1) : 

q a fc-i 



W, 



(2) 



2 5>j 



/>1 



dt 



l+k 



d fe ~ 1 d 2 
( - k + 1) df k + ^ l dt l dt k . l 



(5.9) 



We now express the second additional "ghost-flow" term on the r.h.s. of (|5.7|) as differential 
operator acting on r{t) of a form similar to ( |5.9| ). The starting point are the differential 
equations ( |4.4|) obeyed by the cKP rl r-function, wherefrom we get for the second-order 
derivatives: 

1 d 2 r(t) n ~ l 



T(t) dt r idt rn {=Q 

n-1 l-l 



i=0 j=0 



(5.10) 
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In obtaining relation (|5.10|) we made use of the following Lemma: 
Lemma 5.1 The relation: 

9 -S(f,g) = S(L n (f),g)-S(f,L* n (g)) 



dtr 



n-l 



Y,S(L n - 1 -\$),g)s(f,L* i W 



i=0 



(5.11) 



holds for f an eigenfunction and g an adjoint eigenfunction of the Lax operator L = L rj \ = 
L + + §D~ 1 ^ belonging to the cKP r l hierarchy. 

Proof. We are going to show that d/dt nr S(f, g) = Res (D~ 1 g(L n ) + fD^ 1 ) is equal to the 
right hand side of eq.( |5.11|) (up to a constant). We first apply d/dt mr on the left hand side 
of eq.( |5.11| ). This yields 

d d 



dt dt 



S(f,g) 



d d 



S(f,g) = -R eS (D-\Lr + (9)L m fD- 1 



dt dt 

+ Res(D~ l g L m {L) n + {f)D- 1 
After making the substitutions: 

(L) n + (f) = L n (f)-J2L n - 1 -\<t>)S(f,L* i (^ 



i=0 
n-l 



(L%(g) = L* n (g) + J2L* i ms(L n - 1 -\<!>),g 



(5.12) 

(5.13) 
(5.14) 



i=0 



where use was made of ( |2.21| ), we obtain agreement with the result of applying d/dt mr on 
the right hand side of eq. ( |5.11| ) and using eq. (|2.18| ) as well as Lemma |2j] □ 
Using eqs.( |4.4] ),( |5.10| ) we obtain: 



k-2 



3=0 



J2Uk-2)-j)s(L k - 2 -J(<S>),L*i(y 



k-2 

E 



d 2 T 



2r (*) 1^1 dt r ldt r (k-l- 



•0 



(5.15) 



Collecting ( |5.9| ) and ( |5.15| ), the final form of the cKP r i non-isospectral Virasoro symmetry 
flows reads: 



-X^is-*— 

T 1>1 0t l+r{k-l) 



r(Jfe-l) + l .9 
2r 9t r(fc _i) 



r(fc-l)-l 



+ 7JT E 



d 2 



k-2 



+ ^E 



o 2 



2r dtidt r ( k -i)-i 2 fr{ dt r idt r ( k -i-i) 



r{t) 



(5.16) 



In particular, for cKP^! hierarchies the Virasoro non-isospectral symmetry takes the form 
(eq.flOBD for r = 1) : 



d* k r(t) 



d 



k d 2 

+ E 



^ 1 dti+k-i 2 dt k -i dtidtk-i-i 



it) 



(5.17) 
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Concluding this section it is instructive to point out the relation of ( 5.17 ) with the so 
called Virasoro constraints in conventional discrete matrix models |36], [F| spanning the Borel 
subalgebra of the Virasoro algebra: 



AN) 



o 

oo 

k=l 

oo 



dt 



s > -1 
d 

dt 



s-l 



k+s 



9 a 



dt k dt. 



s > 1 



(5.18) 
(5.19) 

(5.20) 



fc=i 



Here denotes the one-matrix model partition function with iV indicating the size of the 
corresponding (Hermitian) random matrix. On one hand, it can be identified as r-function of 



the semi-infinite one-dimensional Toda lattice model |36]]. On the other hand, from the point 
of view of continuum integrable systems it was shown in |37|, || to be Z N = r^ N '°\ i.e., N-th 
member of the Darboux-Backlund orbit on the subspace of cKPi i r-functions starting from 
the "free" initial r^ 0,0 ) = 1 (see next section for more details about Darboux-Backlund orbits 
on constrained KP hierarchies). Comparing ( |5.19| )-( |5.20| ) with ( 5.17| ) one finds: 



= d* + 



AN) 



d{ + N 2 

C[ N _\ = d* k + (2N + k/2)d/dt k ^ , k>2 



(5.21) 



6 Binary Darboux-Backlund Orbits on cKP rm Hierar- 
chies. Toda Square- Lattice Model 

Let us recall the form of the Darboux-Backlund and adjoint-Darboux-Backlund transforma- 
tions which preserve the constrained form of cKP r m hierarchy Lax operator ( |2.20| ) p7| , |J, 
i.e., we shall discuss anto-Darboux-Backlund transformations for cKP r m hierarchies (for gen- 
eral discussion of DB transformations of generic KP hierarchies without the requirement of 
preserving specific classes of constrained KP hierarchies, see refs. [p8l l8| ) : 

m 

L^L = L + + Y,^D- 1 y* = T a LT- 1 , T a = fc.Dfc- 1 (6.1) 

i=l 
m 

L^l = L + + J2^D- 1 ^ l = f b *- 1 Lf b * , T b = ^ b D% 1 (6.2) 



i=l 



Here $ a , \&& are (adjoint) eigenfunctions entering L__ (p.20 ) with fixed indices a, b which 
henceforth will be assumed a ^ b. Accordingly, for BA functions, the r-function, eigenfunc- 
tions and their respective spectral "densities" , the DB transformations (|6.1| ) imply: 



$ a = T a L($ a ) , *« = ^- 

^ = T a ($ 4 ) , ^ = Tr 1 (^) = -^5($ a ,vl> l ) , i^a (6.3) 
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$BA(t, A) 
<fa{X) 



^T a (^ A (t,A)) 



A 1+ >a(A) 



f (f) = S 8 (f)r(t) 
1 



*.(*) 
^i(A) = A<^(A) 

For the adjoint DB transformations ( |6.2|) we have: 
1 



-A— ^S($ a (t),^A(M)) 

m) = a-v*(a) 



(6.4) 

(6.5) 
i ^ a (6.6) 



$6 

$BA{t,X) 

r BA (t,x) 

^(A) 



* 6 = T 6 L*(* 6 ) 
1 



% ^ b 



-XTr^BA&X)) = X—-rS(^ BA (t,X),%(t)) 



-±T b (r BA (t,X)) ; 
-A 1+ >*(A) ; ipi(X) 



f(t) = %(t)r(t) 
1 



A 



<Pi(\) , #(A) = -Atf(A) 



(6.7) 
(6.8) 

(6.9) 
i ^ b (6.10) 



We shall use the double superscript (n, fe) to indicate the iteration of n successive Darboux- 
Backlund transformations ( |6.3| ) plus k successive adjoint-Darboux-Backlund transformations 
( p.7| ). One can easily show that the result does not depend on the particular order these 
transformations are performed. Therefore, the set of all (n, k) DB transformations, called 
generalized binary DB transformations in what follows, defines a discrete symmetry structure 
on the space of cKP rrn hierarchies corresponding to a two-dimensional lattice. Let us note 



that the (1, 1) binary DB transformations were previously introduced in ref. [2S 
For one-step binary-DB transformed r- and BA functions we get: 



(6.11) 



a s^^^r^t) 

^ (^i°' 0) (^ - [A -1 ]), ^f'°\t - [A -1 ]) 



rBA m (t,x) 



i + 



5($i ' 0) (t),^°' 0) w. 



A 0, (*,A) 



(6.12) 



A 5($l ' 0) (t),^ ' )(t) 

s^Ht + ix-^M^jt + ix-i]) 

s(^°\t)Mr\ t ) 



r BA m (t,x) 



r BA m (t,x) 



(6.13) 



where in the second equalities in (|6.12|) and (|6.13|) we used again (|3.27|) -( [3.28P and Q3.58Q — 
( 3.59 ). Combining eq.( |6.11D with eq.(fOj) for n = 1 we find the following transformation 
formula for the squared eigenfunction potentials: 

d 



(0,0) 



i=i 



8=1 



dt, 



ln(-S(&rKt),¥ b m (t))) (6.14) 
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Let us recall again that here and below a ^ b are fixed indices of Lax (adjoint) eigenf unctions. 
Introducing short-hand notations: 



(6.15) 



we have: 



Proposition 6.1 The following determinant formulae hold for the n-step binary DB trans- 
formed quantities: 



-(n,n 



>(*) 



n-l 



r(°'°)(t) 



?n s *f Jl ( t ).rw=(- 1 ) Bd ? 



$ (n,n) (£) = ( _ 1)f 



s. 



(i-lj-l) 



(6.16) 



det n+ i 


5Sr 1J_1) (*) xr i} (t) 


det n 







(6.17) 



= (_i)«. 



det n+ i 


xl U-D {t) xl (n )[t) 


det n 







s(&r\t)Mr\t) 



det 



n+l 



5. 



(i-lrf'-l) 



det„ 



q(*-1J-1) 



(0 



(6.18) 



(6.19) 



where Xa \ Xb^ an d ^ab^ are defined in ({6.15 ), and the matrix indices i,j run from 1 to n 
or n + 1 according to the indicated sizes of the determinants. 



Formulae (|6.17D — (|6.19 ) can be further generalized to: 



(l m ) s ($i n ' n) (t)) = (-l) n - 



det n+ i 


sV l \t) xt l \t) 


det n 







(6.20) 



n)* 



(n,n) 



(t) 



det n+ i 


xt ( °- l \t) xt {s+n \t) 


det n+ i 


det n 

c(t-i, 

°ab 

q(l+n, 

•Jab 


j - 1] (t) sg+ 


M+n) (t) 
hs+n \t) 


det n 







(6.21) 



(6.22) 
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We are now ready, with the help of (|6.20|) — (|6.22 ), to write down the generalization of 



eq. (|6.16|) for the cKP r-function subject to an arbitrary (n, k) binary DB transformation: 



Proposition 6.2 The general discrete binary Darboux-Backlund orbit on the space of cKP r m r- 
functions, generated by a fixed pair of (adjoint) eigenfunctions $ a , o,nd starting from 
arbitrary "initial" t^ ,0 \ consists of the following elements r^ n ' k ^ : 



det 

fc+i 



-(n,k) 



T 



(0,0) 



S. 



ab 



s. 



(k,j-l) 
ab 



Ao 
Y (*) 



: det 

k 



det 

k+l 



°ab 



-(n-k-1) 



X 



o(<-lj'-l) 
q(n-lj-l) v (n-l) 



An 



for n > k 



(6.23) 



r 



(n,fc) 



(0 



k—n 



det 

n+l 



r (°.°)(t) 



x fe * (n) 



-lrdet 

' n 



s, 



(i-lj-l) 



-(k-n-i) 



X 



det 

n+l 



Q (i-l,k-l) 
D ab °ab 



xt (k-l) 



for n < k 



where vW, v|W and 5^, are as m ( \6.13{ ), and Wk • • • , /*] = det 
caies Wronskian determinant of a set of functions {fi, . . . , /&}. 



1,- 



r? 



(6.24) 



tx,/3=l,...,fc 



indi- 



Remark. Note that the entries in the Wronskians in eqs. fl6.23|) -( p\24T) are determinants 
themselves. 

Remark. In the Appendix we write down the explicit expressions for the cKP rm r-f unctions 
on the most general discrete binary Darboux-Backlund orbit generated via successive (ad- 
joint) DB transformations ( |6.3|) -( |6"7L0|) w.r.t. an arbitrary set of (adjoint) eigenfunctions. 

In the simple case of a "free" initial system with Z/ 0,0 ) = D which, accordingly, is char- 
acterized by: 



-(0,0) 



(t) 



A.a K 'x An 



^Xt ij) = -(-d x ) n xt U) 



(6.25) 



formula fl6.16|) reproduces the Fredholm determinant expression r = det \\5ij + dij\\ for the 
r-function with a^- = f*^ X a i] Xt {i) dy |!J 0. Namely, it follows that: 



d 



1=1 



Res (D-^D^D- 1 



(6.26) 



The latter allows us to identify o^+aj.,- with S (x a l \ xt^^j an d establishes connection between 
the above special case of Q6.16Q and the Fredholm determinant expressions for the r-functions 
of refs.[0, gg]. 

Now, we shall show that the (n, k) binary Darboux-Backlund orbit of cKP r 2 hierarchy 
defines a two-dimensional Toda square-lattice system which describes two coupled ordinary 
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two-dimensional Toda-lattice models corresponding to the horizontal (n, 0) and the vertical 
(0, k) one- dimensional sublattices of the (n, k) binary DB square-lattice. Namely, consider: 



d_ 

dt r 



(n,k) 



(n,k) 



it 11 - ~(n— l,fe) ^(n,fc— 1) 



(6.27) 



where we used ( |6.3| ) and ( |6.7] ). Taking into account the expressions for the (adjoint-)DB 
transformed r-functions ( |6.4| ) and ( |6.9| ), i.e. 



(n,k) 



T (n+l,k) 
q-(n,k) 



(n,k) 



T (n,k+1) 



T 



(n,k) 



eq.(6.27) can be rewritten in the form: 



(6.28) 



Q (n+l,k) (n-l,k) _ (n,fc+l) (n,fc-l) 

ot r 



( T (n,*0) 



or, equivalently: 



(6.29) 



d x — T {n ^ k) - d x T {H ' k) — T ( "' fe) = r («+l.fc) r («-l.fc) _ r (n,fc+l) r (n,fc-l) (6.30) 



Similarly, eq.( |6.27D can be represented as a system of coupled equations of motion for $^ n ' fc - ) 
and vEr^ 1 '*' using again ( |6.28| ) : 



d x — ln$ 

at. 




(6.31) 
(6.32) 



When ty^ 1 '^ vanishes the remaining equations for $^ n,fc - ) reduce for a fixed k to the equations 
of motion for the well-known Toda model on one- dimensional lattice w.r.t. n (and vice versa 
if $("' fc ) = 0). 



7 Discussion and Outlook. Relation to Random Ma- 
trix Models 

In this paper we provided a new version of the eigenfunction formulation of KP hierarchy, 
called squared eigenfunction potential (SEP) method, where the SEP plays a role of a basic 
building block. The principal ingredient of the SEP method is the proof of existence of 
spectral representation for any KP eigenfunction as a spectral integral over the (adjoint) BA 
function with spectral density explicitly given in terms of a SEP. It was pointed out that the 
spectral representations of the (adjoint) BA functions themselves (being particular examples 
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of KP eigenfunctions) can, in turn, serve as defining relations for the whole KP hierarchy 
parallel to Hirota fundamental bilinear identity or Fay identity. 

The SEP method was subsequently employed to solve various issues in integrable hier- 
archies of KP type both of conceptual, as well as more applied character. Many, previously 
unrelated, recent developments in the theory of r-function of the KP hierarchy gained from 
being described by the present formalism. As one of the important illustrations of how 
our method works, we have shown how the SEP, acting on the manifold of wave functions 
i>BA(t,X) by generating non-isospectral symmetry algebra, lifts to a vertex operator acting 
on r-functions. This reproduced in the SEP setting the results of [0, [TJ, [T5|, . 



We have also employed the SEP construction in the context of Hamiltonian reductions 
of KP hierarchy providing: 

• description of the reductions of the general KP hierarchy to the constrained cKP rm hi- 
erarchies entirely in terms of linear constraint equations on the pertinent r-functions; 

• description of constrained cKP rm hierarchies in the language of universal Sato Grass- 
mannian; 

• obtaining the explicit form of the non-isospectral Virasoro symmetry generators acting 
on the cKP^ r-functions. 

The achieved progress should result in further clarification of the status of the cKP rm hi- 
erarchies and their connection to the underlying fermionic field language. It would also be 
interesting to look for the signs of the affine sl(r + m+l) symmetry encountered in construc- 
tion of the cKP r m models by the generalized Drinfeld-Sokolov method associated to affine 
Kac- Moody algebras [|7[ . 

Furthermore, as a principal application, the SEP method was used to derive a series of new 
determinant solutions for the r-functions of (constrained) KP hierarchies which generalize 
the familiar Wronskian (multi-soliton) solutions. These new solutions belong to a new type 
of generalized binary Darboux-Backlund orbits which, in turn, were shown to correspond to 
a novel Toda model on a square lattice. An important task for future study is to find a 
closed Lagrangian description of this new Toda square-lattice model. 

Finally, let us briefly describe another potential physical application of the present ap- 
proach. 



Using the spectral representation for (adjoint) eigenfunctions ( |3.15[) together with ( |2.10|) 



( [2. 11| ), as well as the following form of the Fay identity for r-functions : 



det 



(Aj - Hj)r{t) 



(-1) 



II>j (Aj - Xj) {jm - nj) T {t + *] - HilMi \ 



Tli.j (Ai - fij) r(t) 

(7.1) 

we obtain an equivalent "spectral" representation for T^ n,n ' (t) ( |6.16| ) : 



r M (t) = j n d\ 3 an n - ^) {k - a ■) n (^(a.k^) x 



3=0 i>3 3=0 
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n-1 



j=0 



(7.2) 



Following [41], we can interpret the r-function (|7.2|) as a partition function of certain 
random multi-matrix ensemble with the following joint distribution function of eigenvalues: 



/n—l 
JJ d\j dfij exp {-H(t;{X},{fi})} 
3=0 

H(t; {A}, {/i}) = J2 (Hi(Xj) + #i(^)) + J2 (flaCAi, A,) + F 2 (^, /*,-)) 



+ £# 2 (A i ,// i ) + ff„({A},{/4) 



where the one-, two- and many-body Hamiltonians read, respectively: 

# 2 (A 4 , A,) = - In (A, - Xjf - ln(£ XtX^ 1 ' 3 ) , H 2 (X, //) = ln(A - //) 

^({A}, m) = - \nr^(t + ^[Ar 1 ] - E^r 1 



(7.3) 
(7.4) 

(7.5) 
(7.6) 
(7.7) 



The physical implications of the above new type of joint distribution function (|7.3|)-([T7D 
deserves further study especially regarding critical behavior of correlations. The emerging 
new interesting features of (|7.3|) — ([777D , absent in the joint distribution function derived from 
the conventional two- matrix model ||41| , are as follows: 

(a) the second attractive term in the two-body potential iJ 2 (|7.6|) (both for A- and fi- 
"particles") dominating at very long distances over the customary repulsive first term; 

(b) an additional two-body attractive potential H2 ( [7 -bp between each pair of A- and [i- 
"particles" 

(c) a genuine many-body potential H n ( |7.7|) . 

One of the most important issues here is to exhibit the explicit form of the generalized 
multi-matrix model behind (|7.3|) - (|7.7|) . 



A Appendix: The Most General cKP r m binary Darboux- 
Backlund Orbit. 

Let us first introduce a few convenient compact notations for Wronskian and related Wronskian- 
like determinants: 



W k = W k [fi,...,fk\ =det d a ~ l fp 



W k (f) = W k+1 [h,...j k ,f] 



W k+1 (f) = W[f u ...J k+1 ;f}=det 

fe+i 



s(UJ) S(f k+1 ,f) 



(A.l) 
(A.2) 
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where the matrix indices a, (3 = 1, . . . , k and, as above, d x S (fp, f) = fp f ■ The Wronskian- 
(-like) determinants ( |A.l|) -( |A~2j ) obey the following useful identities: 







'Wk-i(/A w k (f)w, 



k-l 



w k 







' w k+1 (fy 



W k (f)W, 



k+1 



(A.3) 



where the first one is known as Jacobi expansion theorem (see, e.g. |fl2|), whereas the second 
identity in ( A.3| ) can be easily verified via induction. Eqs.( A.3) imply in turn the identities: 



T, 



with: 



j 



Ti(f) = 



w. 



*-l 



D- 



rri* — 1 
3 



Tr\f) 



Wkjf) 



(A.4) 
(A.5) 



Now we can use eqs. ([A.4|) to derive explicit expressions for the (adjoint) eigenfunctions 
and r-functions of cKP rjm hierarchies generated via successive (adjoint) DB transformations 
( p.3| ) — ( |67TO| ) w.r.t. an arbitrary set of (adjoint) eigenfunctions of the "initial" cKP rm Lax 
operator L = L rm Q2.20|) . We shall denote the latter arbitrary successive (adjoint) DB 
transformations by the following double- vector superscript: 



(ft, kj = ((m, . . . , n m ), (kx, . . . k m )) 



(A.6) 



indicating ri\ successive DB transformations w.r.t. $i etc., until n m DB transformations 
w.r.t. $ m and, similarly, k\ successive adjoint-DB transformations w.r.t. ^ 1 etc., until k m 
adjoint-DB transformations w.r.t. \P m . Specifically, we have: 
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(o) 
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(ni-l) 
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y (0) 
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X 



(n m -l) 



n,0 
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r (o) 
Xi , 
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,Xi 
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v (0) v («m- 

Am ) ■ ■ ■ ) Am 




(0) 

Xi r 


■,xl ni 


-1) 


y (0) 




!-2) 


JO) (nm-l)' 


\V 


(0) 




-1) 


Y (0) 


v ( ™ 

■;Aa 




.,(0) v (nm-l) 



, for n a > 1 



w 


" (0) (0) (n m -l).^(5,B)" 


w 

V 


V (Q) yK" 1 ) v (0) v ("m-l) 

Al ) * • ")A,1 • 5 Am 5 ■ ■ ■ ) Am 





(A.7) 
(A.8) 

(A.9) 



for n„ = 



where the functions xffl are the same as in (|6.15|) with the superscripts (0, 0) replaced with 
the corresponding double- vector ones (o, Oj . Eqs. (A.7)-(Q?) already appeared in |28| (see 
also refs.|f43f), whereas eq. (|A.9|) is derived via iterative application of the second identity in 
( |A.4| ) and taking into account ( |A.7| ). 

Now, performing arbitrary successive adjoint-DB transformations on r( n '°) (|A.8|) accord- 
ing to the second eq. ( |6.9D upon using the first identity in ( |A.4j ) and inserting there the explicit 
expressions ( |A.9f) , we arrive at the following: 
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Proposition A.l The most general discrete binary Darboux-Backlund orbit on the space of 
cKP^ T-junctions is built-up of the following elements: 



(0,0) 



-W 



v (°) v (rtl_1) v(°) v^" 1 ) 

Al j • • • ) Al 5 • • • 5 Aa ! ' ' ' ' Aa 



" Ea+l k i 



W 



An An An An 

^(0,o+l)) • • • ) ^(fco+i-l.o+l)' ' ' ' ' ^(0,m)) • • • ) ^(fem-ljm) 



(A.10) 



where: 



v (°) v (ni_1) v(°) v K-l)-v* {0 

Al j • • • ) Al ) ■ ■ ■ j Aa ' ' ' ' ' Aa ' As 



(n, £) = ((ni, . . . , n a , 0, . . . , 0), (0, . . . , 0, k a+1 , ...,k m )) 
and, furthermore, notations (6.1d[) and (A. 2) are employed. 



a = 0, 1, 



, m 



(A.11) 



(A.12) 



Remark. The reason for the zero entries in the labels (|A12|) of the most general binary DB 
transformations, preserving the spaces of cKP r?n hierarchies (|2.20 ), lies in the fact that any 
pair of two successive (adjoint-)DB transformations w.r.t. $ a , \I/ a , i.e. both with the same 
index, is equivalent to an identity transformation as one can easily conclude by combining 
the second equation in ( |Q|) with the second equations in ( |6.4| ) and (|6] 
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